Hyperon ordering in neutron star matter 
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Abstract. We explore the possible formation of ordered phases in neutron star matter. In the 
£NJ ■ framework of a quantum hadrodynamics model where neutrons, protons and Lambda hyperons 

interact via the exchange of mesons, we compare the energy of the usually assumed uniform, liquid 
phase, to that of a configuration in which di-lambda pairs immersed in an uniform nucleon fluid are 
■ localized on the nodes of a regular lattice. The confining potential is calculated self-consistently as 

resulting from the combined action of the nucleon fluid and the other hyperons, under the condition 
of beta equilibrium. We are able to obtain stable ordered phases for some reasonable sets of values 
of the model parameters. This could have important consequences on the structure and cooling of 
neutron stars. 
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^ ; The equation of state in the interior of neutron stars is usually considered to be that 

of an interacting Fermi liquid in a uniform phase. Under certain conditions however, a 
crystallized phase may be energetically more favorable. This in turn can have interesting 
consequences on the structure and evolution of the neutron stars. We may quote as 
examples thereof (i) triaxial configurations with emission of gravitational waves []Ip, 
(ii) modification of the oscillation modes of the star or (Hi) of the neutrino transport 
properties 0. 

In a previous paper [Q] we presented the results pertaining to a simplified model 
where we considered three baryonic species: the neutron, proton and Lambda hyperon, 
interacting with each other through the exchange of a and £0 mesons. We found that 
it was possible to find some sets of the model parameters such that the energy of a 
configuration where pairs of A hyperons with antiparallel spins located at the nodes of a 
cubic lattice, could be energetically more favorable than the corresponding liquid phase. 

The model presented in [Q] was a strictly minimal one. Here we would like to perform 
again these calculations with a more realistic description of the nuclear interaction. In 
order to reproduce the correct values of the incompressibility modulus and the asymme- 
try energy, we choose an approach based on the Density Dependent Hadron Field theory 
(DDRH) developped by the Giessen group [@, 0]. In this model, the nucleon-meson 
vertices have a functional dependence on the density operator. In this way, some basic 
features of the short-range correlations are taken into account. 



MODEL OF BARYONIC MATTER 



Lagrangian 

We describe the nuclear interaction in the framework of a phenomenological quan- 
tum relativistic model. The interaction piece of the Lagrangian density in the Density 
Dependent Hadron Field Theory (DDRH) reads: 
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where the couplings Y aB are density-dependent through their functional dependence on 
the baryonic current . 
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We will consider only one species of hyperons, the A. In a realistic model, it can be 
expected that the E hyperon will also play a role. However, the interaction of this hyperon 
with nuclear matter is less well known, so that it has not yet been implemented in the 
framework of a DDRH-type model. 



Liquid phase 

In the liquid phase we impose the usual conditions for charge neutrality and chemical 
equilibrium under /3 decay among the baryons (neutron, proton, A hyperon) and leptons 
(electron, muon) taken into account in our model. 

n e J r n il =n p , \l n = \l K , \L p - \L n = fl = \L e = 

These equations determine the relative fraction of each species of particles. The equation 
of state is then obtained from the energy momentum tensor. For details we refer the 
reader to [@, 0]. I n particular we will be interested in comparing the energy density of 
the liquid phase p L = T 00 to its counterpart in the crystallized phase. 

Self consistent confining potential 

We now assume that the hyperons are ordered on a lattice and are localized in gaussian 
clouds. We further assume that there are two hyperons with antiparallel spins per lattice 
site . This is partly in order to avoid complexities related to spin-spin interaction, but we 
can also justify this choice by arguing that the A-A interaction, as inferred from data on 
double hypernuclei, is attractive and may even favor a bound state (~ H-dibaryon). In the 
spirit of the Sommerfeld aproximation, we neglect the redistribution of the surrounding 
nucleons. 



The potential energy of an hyperon around a lattice site is calculated as the sum 
of the interaction potential energies generated by hyperons at other lattice sites r- = 

a(li + mj + nk) . 

t/su P (r)=£[/(r-F.) (3) 
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In this expression, the potential energy of each hyperon is obtained by taking the 
convolution 

171 f\ r ~\~R\ 

U(T) = — Rn(R)dR V 0BE (x)xdx (4) 
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of a gaussian distribution 

-^°J e- M ^f 2 (5) 

with the elementary one boson exchange potential V 0BE (x) corrected for the finite size 
of hyperons by form factors (with cutoffs A a , A a ) and for relativistic effects in the 
momentum expansion. Note that in the present configuration, only the central and spin- 
spin components contribute. 

The potential calculated in this way is then approximated by a parabola £/ par parame- 
terized by its depth U Q and the frequency of oscillation of the hyperons in the potential 

well v = yJV 2 U pai (T*)\y>_ Q /M A . The width T = l/y/M A v Q of the gaussian distribu- 
tion @ is determined self-consistently by the frequency of this harmonic oscillator. 



Equation of state of the ordered phase 

We again impose the conditions for /3 equilibrium 

\l n = \l p + \l e , [l n = fl A , [l e = [1^ (6) 
with the chemical potentials given by 
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The nucleons evolve in constant mean fields consisting of two contributions: To the field 
produced by the homogeneous nucleon fluid (e.g. ft) ), we add the spatial average of 
the fields generated by the periodic hyperon distribution (e.g. < ft) A > s ). The chemical 
potential of the hyperons is now determined by the parameters of the confining potential. 
The rearrangement terms E rearr contain derivatives of the couplings Y aB with respect to 
the density. 



We further impose the conditions of charge neutrality +n e = n p , and express the 
baryonic density n B = n n + n p + n A as the sum of the nucleon densities and the hyperon 
density, the latter being related to the lattice parameter a by n A = 2/a 3 . Finally, we write 
the defining equations for the effective baryon masses and the self consistent equations 
obeyed by the fields o, ft) , p° and 8. 
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The equations for chemical equilibrium again determine the chemical com- 

position of the matter at a given baryonic density. They have to be solved self consis- 
tently with the equations for the confining potential to which they are interrelated 
through the effective hyperon mass, M A the lattice parameter a, the depth U Q and the 
oscillation frequency v . 



NUMERICAL RESULTS 



We will use the rational parametrization of the coupling constants suggested in Ref. 
with the parameters of Ref. [[J] 
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(9) 



This function is chosen so as reproduce parameter free Dirac-Brueckner calculations in 
the nucleon sector. The parameters x aB are adjusted in order to reproduce the binding 
energy of the A in hypernuclei. The coupling constants used in the present work corre- 
spond to the "model 1 with DD phenomenological" parametrization of Reference [[7p. 
We would like to stress the fact that, in contrast to the calculation performed in [Q], we 
have essentially no free parameters in the results presented in this work. 

When numerical convergence is reached, we obtain the parameters U , v , M ( , p Fi , 
etc. which fully determine the thermodynamical state of the system. We now are in a 
position to calculate the energy density of the crystallized phase. It is obtained as 



P C = P P +Pn + Pe + Pfi+P A + P fields 



(10) 
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Figure 1. Energy difference between the 
liquid and crystallized phase, as a func- 
tion of baryonic density 



Figure 2. Pressure in the liquid and the 
crystallized phases, as a function of bary- 
onic density. 
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This energy is then compared to the energy density of the corresponding liquid phase 
p L . The result is displayed in Figure 1. With the chosen model parameters we obtained 
an ordered phase which is energetically more favorable than the liquid one above the 
threshold for hyperon production n lh = 1 .95 n sat . The shape of this curve is similar to the 
one which was obtained for parameter set C in with however a somewhat smaller 
energy gain. The stability of the solid phase increases at high density. 
Finally, we obtain the pressure from the thermodynamical relation 

P = i#*) (ID 
dn B 

The result is displayed in Figure 2. In contrast with the results obtained for parameter 
set C in [f|], the equation of state of the ordered phase is now slighlty harder than the 
liquid one. 



CONCLUSION 

We investigated the possible formation of an ordered phase of the baryonic matter 
present in the core of the neutron stars, in which the hyperons are located on the nodes 
of a cubic lattice. The model presented in this short contribution is a first upgrade 



of our original barebones model [|]]. We now have a more realistic description of 
the nuclear interaction, in particular of the nucleon sector. Our main conclusion is 
that the findings of [Q] are confirmed, namely the ordered phase can be energetically 
more favorable for some choices of the model parameters. As could be expected, our 
result is found to depend more strongly on the features of the - largely unknown - 
in-medium hyperon-hyperon potential than on the properties of the nucleon background 
such as incompressibility modulus and asymmetry. A more systematic study with several 
alternative parametrizations is underway and will be presented in a further publication. 

The model which we have presented here can be improved in several ways. In order to 
be consistent with the DDRH picture, we described the hyperon-hyperon interaction by 
sigma and omega exchange with density dependent couplings in the form T Jn N , n») = 
x aA T aN (n N + n A ) with a G {a, Co}. On the other hand, some authors also introduce 
additional o* and exchange in order to reproduce the attractive A — A interaction 
inferred from the available data on double hypernuclei. Furthermore, parametrizations 
of the free hyperon-nucleon and hyperon-hyperon potentials consider the exchange of 
e.g. kaons or r\ which are also not present in our calculations. 

Ideally one could extract an in-medium hyperon-hyperon potential from Brueckner- 
type calculations and test the result comparing with the available hypernuclei data. 
Whereas several non relativistic calculations exist in the litterature [ [T0| , |TT| ], the corre- 
sponding calculations have not been performed in a relativistic framework. In particular 
a full fledged relativistic calculation including all hyperons (not only A but also E and S) 
and a subsequent DDRH-type parametrization would be especially welcome. It would 
allow us to have a greater variety of lattice configurations involving in particular the E~ 
hyperon. Indeed the E~ is formed in the liquid phase at the same density as the A in 
most neutron star matter calculations. 

An other direction of research would be to release the Sommerfeld approximation 
and take into account the redistribution of nucleons. This involves considering screening 
correlations of the RPA type to the potential. 
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